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Abstract 

N— level quantum model is proposed in which the energies are represented by 
an N— plet of zeros of a suitable classical orthogonal polynomial. The family 
of Gegenbauer polynomials G(n, a, x) is selected for illustrative purposes. 
The key novelty lies in the use of non-Hermitian (a.k.a. cryptohermitan) 
Hamiltonians H ^ W . This enables us to (1) start from elementary secular 
equation G(N, a, E n ) = 0, (2) keep our H, in the nearest-neighbor-interaction 
spirit, tridiagonal, (3) render it Hermitian in an ad hoc, non-unique Hilbert 
space endowed with metric 8 ^ /, (4) construct eligible metrics in closed 
forms ordered by increasing nondiagonality and (5) interpret the model as a 
smeared iV— site lattice. 



1 Introduction 



In atomic, molecular, nuclear and solid-state physics the simulation of quan- 
tum phenomena via finite- dimensional Schrodinger equations 
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is often motivated numerically. Indeed, whenever a realistic Hamiltonian gets 
approximated by its suitable N by N simplification H = H^ N \ the numerical 
solution of Eq. ([1]) becomes routine p] , especially when our finite-dimensional 
Hamiltonian is chosen tridiagonal, 
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At a fixed iV, various N by N matrix problems (CQ) + ([2]) are often used in 
ambitious phenomenological considerations since matrices (sometimes 
reinterpreted as the so called chain-model Hamiltonians or lattice Hamil- 
tonians with the nearest-neighbor interaction) may mimic, say, a solid-state 
spectrum of energies in condensed-matter physics. These concepts found new 
applications in the context of apparently non-Hermitian versions (we would 
rather call them "hiddenly Hermitian" or "cryptohermitian" [2] versions) of 
the XXZ spin chains [3], of the Bose- Hubbard models [4], of the Friedrichs- 
Fano- Anderson tight-binding lattice models [5] , of the tightly bound lattices 
of electrons [6], optical lattices [7] etc. There exist many other papers which 
are also certainly worth mentioning. In their incomplete complementary sam- 
ple we would like to attract attention of the reader to the close connections 
between non-Hermitian chain models and the so called Reggeon field theory 
[8] or to Ising model and quantum spin chains [9]. 



The combined mathematical and physical appeal of the generic discrete 
and tridiagonal models (j2j) seems partially marred by the more or less purely 
numerical [TU] or perturbative [TTJ nature of their solution. For this reason, 
analytically solvable models are often preferred in analysis [12] . Some authors 
simplified mathematics by paying attention to the effects connected with 
the restricted, one-parametric variation of the end-site-interaction matrix 
elements Cj and bj + i with j — and j = N — 2 in Eq. (j2]) [6j [13] . 

Marginally we could add that similar discrete solvable models with pairs 
of point-like interactions played important role in the recent extensive dis- 
cussion of some conceptual problems of cryptohermitian quantum scattering 
[TU IT5] . In this context a lot of misunderstandings emerged when people for- 
got to distinguish between the "formal coordinate" x (often chosen as playing 
the role of the argument in wave functions ip {%) ) and the "observable coordi- 
nate" (a position-operator eigenvalue denoted by another symbol, say, q). In 
a very well written paper [16] interested reader may find the nice explanation 
of this subtlety emerging as highly relevant even on the very elementary level 
of mathematics used in introductory textbooks on quantum mechanics. 

Once one moves to the more sophisticated cryptohermitian models where 
the "formal coordinate" x itself ceases to be observable, the concept of "lo- 
cality" must be reconsidered and used with enhanced care. For example, a 
very instructive comment given in section 5 of Ref. [17] shows that the formal 
wave function of a physical localized state may look non-local as a function 
ip{x) of the formal coordinate x. 

In opposite direction it has been noticed and emphasized by Jones [H] 
that in virtually any experimentally oriented setup we usually treat interac- 
tion V as if it were prepared as a specific function of the measurable coordi- 
nate q. In this sense, the crucial role of the specification of observables and 
of the difference between x and q gets even more important in non-Hermitian 
setting [U]. 

In order to circumvent similar complications a number of papers studied 
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just bound-state problems and preferred their exactly solvable non- Hermit ian 
models [19J. The solvability-guaranteeing simplifications may reduce the 
menu of interesting phenomena. Typically, the simplified models explain 
the emergence of fragile, unstable components in the spectra [20] but they 
can hardly compete with realistic models in offering sufficient variability of 
the parametric dependence of the energies [3]. The spectra obtained in the 
simplified solvable model of Ref. [6] admit, for example, just a very special 
form of the confluence of energy pairs while a much richer menu of quantum 
catastrophes of this category may exist in general [2Tj . 

A remedy has been found in Ref. [22]. We revealed that there exist non- 
numerical chain models or quantum lattices (j2J) with a much less restricted 
qualitative variability of spectra. These models were characterized by a de- 
localized interaction exhibiting an up-down symmetry. The pairs of sites 
with indices m and N — m were attached the same strength of impurity or 
interaction. Although the productivity of such an artificial assumption was 
reconfirmed, say, in refs. [B] and [23], its physical interpretation remained ob- 
scure. One would like to have some exactly solvable quantum-lattice models 
without such a symmetry. This motivated our present analysis during which 
we developed another class of solvable quantum-lattice models of form (T2]) 
without similar non-local, long-range auxiliary correlation. 

2 Gegenbauer-polynomial quantum lattice 

In connection with the definition of the concept of solvability misunderstand- 
ings frequently emerge. The puzzle may find different resolutions. In a 
context-dependent way the property of being solvable is assigned, e.g., to 
differential Hamiltonians H = p 2 + V(x) for which all of the wave functions 
(x\^ n ) of bound states prove proportional to suitable classical orthogonal 
polynomials [23]. In our present paper we shall transfer such a definition of 
exact solvability to the difference and finite-matrix equations. Thus, we shall 
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postulate that the iV— plet of our iV— dimensional bound-state vectors \i/j n ) 
in Eq. ([T]) is given in advance. 

Naturally, the most straightforward definition of these vectors would spec- 
ify them directly in terms of some classical orthogonal polynomials. For 
the sake of brevity we shall solely pay attention to Gegenbauer polynomials 
G(n,a,x) (= C*(x) in [25] or Cn\x) in [22]; our notation is taken from 
MAPLE [27]). As long as these (sometimes called ultraspherical) polynomi- 
als degenerate to the different (viz., Chebyshev) polynomials at a = 0, we 
shall assume that a > 0. In this case they satisfy the well known recurrence 
relations 

n G(n, a,x) = 2 (n + a — 1) x G(n — 1, a, x) — (n + 2 a — 2) G(n — 2, a, x) (3) 

at n — 1, 2, . . ., with initial G(0, a, x) — 1 and G(l, a,x) =2 ax. 

In the initial step of our constructive considerations we shall guarantee 
the validity of our above-mentioned matrix Schrodinger Eq. ([T]) by assuming 
its formal coincidence with the truncated version of recurrences (J3]). This 
means that we shall just use the following input form of the bound-state 
eigenvector, 



I (0\^ n N) ) = G(0,a,E n ) \ 
(l\^ n N) ) = G(l,a,E n ) 



\ (N-l\4 N) )=G(N-l,a,E n ) J 

and determine the n— th energy level E n as the value of coordinate x at which 
recurrences ()3]) terminate. Thus, every energy will coincide with one of the 
roots of the closed-form secular equation 

G(N,a,E n ) = 0. (5) 

Our Gegenbauerian Hamiltonian H = H(a) will just mimic recurrences 
Its main diagonal will vanish (i.e., we set ao = a\ — . . . — in and the 
pair of non-vanishing neighboring diagonals will be composed of elements 
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numbered by j — 0, 1, . . . , N — 2, 

c . = Cj ( a ) = l/(2a + 2j) , b j+1 = b j+l {a) = (2a + j)/(2a + 2j + 2) . (6) 

This idea forms the starting point of our abstract message in its concrete 
Gegenbauer-polynomial realization. Within the more general class of quan- 
tum lattices and discrete models ($Z§ exemplified by such a choice the matrix 
elements are real but the matrix H itself is, generically, asymmetric, i.e., 
non-Hermitian. Fortunately, its spectrum is real (i.e., potentially observ- 
able) so that we are allowed to treat this H as an exactly solvable effective 
Hamiltonian of a quantum system with the prescribed segment of spectrum 
fitted by an N— plet Ei t N \a) of roots of Gegenbauer polynomial G(N, a, E). 

3 (Hidden) Hermiticity 

It is known that the manifest non- Hermiticity feature does not disqualify 
operator H ^ W from being used as a Hamiltonian of a quantum system. 
After all, not too dissimilar non-Hermitian phenomenological Hamiltonians 
(complex and acting in a finite-dimensional vector space) were used in Refs. 
[5] - [7]. Interested reader may find a compact introduction into quantum 
theory with similar cryptohermitian Hamiltonians either in our review [2] or 
in this section. 

In essence, we must get rid of the overrestrictive and most elementary 
(often called "Dirac's" [28J) requirement of the current but very special Her- 
miticity defined via the mere vector or matrix transposition accompanied by 
complex conjugation. This defines dual vectors called, in the conventional 
textbook language, "Dirac's bra- vectors" , 

r (Dirac). |^ ^ (7) 

The choice of , ]~( Dtrac ) (represented just by appended superscript t when 
applied to operators) is not the only option. In models with Dirac-non- 
Hermiticity H ^ we must necessarily use another, less trivial definition 
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of Hermitian conjugation. The point is that after such a change of definition 
our operator H may become self-adjoint and compatible with postulates of 
Quantum Mechanics. 

The transition to general Hermitian conjugation will require a modifi- 
cation of conventional notation. Firstly, the "new" dual vectors must be 
defined by generalized formula 

r (e) : m -> «v>h=Me (8) 

where matrix G is called "metric" (18] and where, whenever ^ J, the 
resulting dual vectors are marked as "brabras". Secondly, the same danger 
of misunderstanding threatens the application of the non-Dirac Hermitian 
conjugation to operators A so that we recommend it to be marked by a 
different (viz., doubled) superscript, 

A a % := e- x A f e. (9) 

In the spirit of any good textbook on Linear Algebra, Functional Analysis 
or Quantum Mechanics the metric must be required invertible, Hermitian 
and positive definite [18] . After two notation innovations jSJ) and (Q the 
formalism of Quantum Theory remains unchanged. On the level of nota- 
tion the symbol of double bras ((() will replace all the Dirac's simple bras 
((), especially whenever a mean value, physical probability or measurements 
are concerned. Similarly, in formulae carrying physical meaning the simple 
superscripts ' must be all replaced by their doubled forms The "false" rep- 
resentation H( F > of the Hilbert space with the Dirac's unacceptable 0^ = I 
must consequently be replaced by the "standard" Hilbert space of phys- 
ical states if). 



4 Hilbert-space metrics 



From the pragmatic point of view the theoretical imperatives of preceding 
section may be softened, during practical calculations, by staying in the naive 
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(and, by assumption, much friendlier) Hilbert space and by the treat- 
ment of the obligatory doubled bras (( of Eq. (jSJ) and doubled superscripts * 
of Eq. OH]) as mere abbreviations. In Ref. [2] we summarized further reasons 
for a parallel use of spaces together with their "friendly - false" partners 
H( F ) . Firstly, just the knowledge of the matrix (which must be self-adjoint 
in [18J) is fully sufficient for all purposes. Secondly, the key guarantee 
of unitarity of the evolution generated by H in T-L^ (where H = H$ and 
8^J) gets easily legible in via "translations" (JSJ) and flUJ). In fact, a 
deeper explanation of this point deserves the (following) separate paragraph. 

4.1 Dieudonne equation 

In Ref. [2] we explained the way in which the Hermiticity of H in (based 
on the nontriviality of metric ^ I) may be understood as equivalent to 
the manifest Hermiticity of a suitable isospectral operator 

f) = VlHn- 1 = f) f . (10) 

The latter operator is defined, in principle, in another, third Hilbert space 
*HS p ' with trivial metric = I (the superscript stands for "paternal" or 
"physical"). It is assumed that spaces and are unitary equivalent 
so that we may recall Eq. ([9]), deduce 

f) 1 " = (fi- 1 ) t #tQt ; (n) 

abbreviate Q^Q := G and end up with the relation 

H^Q = QH (12) 

dating back to the old paper by Dieudonne [29J. That's why we shall call 
Eq. (I12p "Dieudonne's equation" in what follows, keeping in mind that this 
is meant in a loose sense since Dieudonne himself admitted that in (TP2]) 
might not be invertible. 

For our finite-dimensional real Hamiltonians H = H^ N ' which are given in 
advance, the latter equation forms the set of N 2 constraints imposed upon the 
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[N(N + l)/2]— plet of the unknown real matrix elements of matrix = 0t. 
Our task may now be formulated as a non-numerical construction of complete 
solution of this linear algebraic system. 

4.2 The method of solution 

The constructive way of making Hamiltonian H and metric compatible 
with Dieudonne's Eq. (fl2|) is not too easy in general. The main result of 
our paper will be the non-numerical construction of the general metric 
which will satisfy Eq. ( I12p for the Gegenbauerian input Hamiltonian H^ N \a). 
Ipso facto, this will also make our Hamiltonian self-adjoint in the respective 
physical Hilbert space T-L^ s \ 

In full detail, the construction of metrics will be described in section 
below. In a preparatory phase let us now just explain its key ideas. Firstly, 
in the light of the linearity of Eq. (Tl2|) we shall assume that the metric may 
be sought in the form of superposition of certain simpler matrices V which 
will satisfy the same equation, 

{H^ N \a)) ] V = VHW(a), (13) 

but which will not necessarily be invertible or positive definite. Secondly, 
we shall assume that these "pseudometric" matrices will form an iV— plet of 
linearly independent solutions V = v\. (a) with k = 0, 1, . . . , N — 1. This 
will enable us to search for the metric in the form 

iV-l 

= eW( a » = oqQW ( a ) + akVl N) (a) (14) 

fc=i 

where the variability of the N— plet of real parameters a = (ao, ati, ■ ■ ■ , ckjv-i) 
will only be restricted by the requirement of the positivity of the matrix 

The concrete implementation of the requirement of the simplicity of the 
individual auxiliary pseudometrics Vk is model-dependent. For our present 
model their explicit construction proved feasible when we assumed that every 
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Vk is a (2k + 1)— diagonal matrix. This assumption itself resulted from the 
experience which we gained during the similar constructions of metrics as 
performed in Ref. [30] • This experience also facilitated the organization of 
our concrete recurrent calculations. 

The key idea of our present non-numerical algorithm of solution of Eq. ( Fl3l) 
remained the same as in Ref. [3U]- In concrete applications we shall see 
how this recipe employs the chess-board-like "coloring" of elements of rele- 
vant matrices. In this manner, each Hamiltonian H gets separated into its 
"white- field matrix elements" (say, all elements Hj^ with \j — k\ =even) and 
"black- field matrix elements" (i.e., elements Hj^ with \j — k\ =odd). Once 
the same coloring is applied to the ansatz for the metric (or rather to 
each indefinite and sparse pseudometric ( Q ))> one is immediately able to 
decompose Eqs. (TT^|) and/or (TT5j) into their "same-color" subsystems and to 
develop and employ some suitable ansatzs for their recurrent solution. 

One should not forget that even before finishing the systematic construc- 
tion of all of the components v\, (a) of the metric we may interrupt the 
process and turn attention to the truncated versions of series (|14|) . 

e£V, a ) = + E «;*f°(«) • ( 15 ) 

3=1 

Here, the mere k free parameters <x, appear arranged in a shorter, primed ar- 
ray a'. One should also pay attention to the fact that in Ref. [3U] as well as in 
our present model the k— subscripted special metrics (1151) still remain sparse, 
containing just 2k + 1 no n- vanishing diagonals. The latter observation will 
certainly facilitate our ultimate task of imposing the positivity requirements 
upon expansions (|14p or ffTBl of the metric. 
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5 Diagonal metrics 



5.1 The construction of Go (a) 



All the details of the implementation of our above recipe depend on the form 
of the input Hamiltonian H. For its Gegenbauerian choice given by Eq. (jHj), 
this Hamiltonian is an extremely elementary, purely "black-field" matrix, 
rendering the recurrent solution of Eq. ( fl2|) particularly straightforward. For 
illustration purposes let us now consider the diagonal (i.e., k = 0) ansatz 



Oo 

0i 

... 














9 N _ 2 

. . . o e N -i 



(16) 



As long as the individual matrix elements will not vary with the growth of 
dimension N, we may leave the value of N unspecified. The inspection of re- 
currences (|12p then reveals that they connect just equal-color elements. This 
means that a priori, ansatz f[T5"j) may lead to nontrivial solutions. We may 
start their recurrent construction from any nonvanishing element, say, from 
9q = 2a 2 . After a comparatively tedious algebra this choice of normalization 
leads to the compact and transparent final result with 9i = a + 1 and with 

a + j 



0; 



;i + 2a)(2 + 2a)...(j-l + 2a) 



(17) 



at all the remaining j 



2,3,..., N- 1. 



5.2 A comment on matrices f) 

An important feature of the above-constructed metric Go (a) is that it is easily 
invertible and manifestly positive definite at any a > and at any N > 1. 
The existence of such a metric is an important merit of the model because 
we may now recall relation f lTU]) . define the matrix elements of the simplest 
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auxiliary matrix = Qq, 

(^o) mn = 6 mn (18) 
and obtain finally the simplest explicit partner Hamiltonian 
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acting in space = Hq , isospectral with our original non-Hermitian 

matrix H^ N \a) and possessing matrix elements easily derived in closed form, 

1 / 2a 4- k 

^2\j {a + k)(a + k + l y * = o,i,...,w-2. (20) 

Due to the unitary equivalence between Hilbert spaces U { P and = H { S) 
we may conclude that Eq. ( !T6|) represents the simplest possible Hermitization 
of our Gegenbauer-oscillator Hamiltonian H^ N \ 

The existence of the partner Hamiltonian f )19p trivially re-confirms the 
well known fact that the spectra of energies E { n \a) defined by Eq. © arc 
all real [26]. Moreover, the manifest positivity and diagonality of ®o{a) makes 
the explicit construction of matrix Qq virtually trivial. The latter observation 
is not easily transferred to other models. For example, interested readers may 
consult Ref. [17] showing that and in which way a very simple Hamiltonian 
H may be assigned extremely complicated isospectral partners f). 

Exceptions from the latter generic rule exist. In the present context of 
models on lattices a typical one has been found in paper [21]. A non-diagonal, 
band-matrix metrics O has been shown there to admit a transparent, sparse- 
matrix structure of factors in 6 = Q^Q as well as of the corresponding 
isospectral Hamiltonian f). Of course, this type of result must be considered 
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exceptional. Formally, the reason is that the use of formula f JTOj) which defines 
the partner Hamiltonian f) requires the explicit knowledge of the inverse 
matrix Q~ l which is usually not a sparse matrix even if Q itself is. 

This being said it is necessary to admit that one cannot exclude that 
our present Gegenbauerian example will prove exceptional and that it will 
also admit the existence of compact formulae for f), e.g., at some coordinate- 
smearing choice of k — k( except ionai) > 1- With the notable exception of our 
knowledge of tridiagonal k( except i ona i) = matrix ffT9l the existence and possi- 
ble structure of such formulae is an open problem at present. In fact, the lack 
of our explicit knowledge of all of the manifestly Hermitian Hamiltonians f) 
hinders, first of all, the most common strategy of interpretation of the sys- 
tem in question illustrated, e.g., in Ref. [10] and based on the correspondence 
principle applied directly inside / W- P >. 

A positive aspect of the existence of missing parts of the puzzle is that 
if any relevant matrix f) really remained sufficiently simple and defined in 
closed and compact form, all the reasons for working with its equivalent 
representation H in would in fact be lost. The situation is similar to the 
preference of H in nuclear physics [18] where the more complicated partner 
f) is even well known in advance. The same preference of the maximally 
simple representation of the Hamiltonian remains recommended for concrete 
calculations even though we proceed here in opposite direction, viz, from 
the choice of H to the construction of its Hermitizations mediated by 0s in 
alternative Hilbert spaces 

6 Band-matrix metrics 

It has been explained in Refs. [15] and [30] that tridiagonal metrics, i.e., in 
our case, the one-parametric family of matrices 

e[ N) ( ai ,a) = Q { N \a) + ai v[ N \a) (21) 
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simulate a nearest-neighbor smearing of coordinates while the pentadiagonal 
metrics 

iW^, n . n uoWux n ,.#)/ n u n ,-DWu (22) 



e^ ; (a ls a 2 ,a) = 6^(a) + Ql 7^ » + a 2 Vf'(a) 



may mimic a next-to-nearest neighborhood smearing, etc. In this manner 
the index k in Eq. ffl5|) is tractable as a certain measure of a dynamical, 
Hilbert-space-related "nonlocality" of the quantized lattices in question. 



6.1 Tridiagonal metrics 0[ N \ai,a) 

In Gegenbauer example ([6]) all the generalized k = 1,2, . . . metrics (fl4|) may 
be constructed in closed form, non-numerically, by the recurrent solution 
of Eq. f )12p . After some trial-and-error experimenting the first nontrivial, 
tridiagonal metric 0j N \a\, a) (containing just the single item in the primed 
array of parameters 3 = cti) may be found via the tridiagonal (or, more 
strictly speaking, bidiagonal) k = 1 ansatz for its only nontrivial sparse- 
matrix component 
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(23) 



The combined use of the experience and computer algebra leads to the 
truncation-independent result. Using the convenient initial Ki = 2a and 
k,2 = 1 one obtains the closed formula 

1 



(24) 



J (l + 2a)(2 + 2a)...(;'-2 + 2o) 

for the solution flU} of Eq. (JT2J valid at all j = 3,4, . . .,N - 1. Let us 
re-emphasize that these matrix elements exhibit the remarkable property of 
not changing their form with the matrix dimension N. 
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6.2 The domains of positivity of metrics 0| \a>i,a) 



It is worth noticing that the positive definiteness of the tridiagonal metrics 
(l2Tj) would be lost for larger a\ > oi^ tical (a). Using an analytic method this 
expectation may be illustrated via a slightly renormalized two-dimensional 
metric 

2 a 2 ab 



6< 2) (&/2,a) 



possessing two real eigenvalues 



ab a + 1 



1/2 a + 1/2 + a 2 ± 1/2 ^ -3 a 2 + 2a-4a 3 + l+ 4a 4 + 4 (ab) 2 . 

It is easy to deduce that the domain D of positivity of this metric coincides 
with the interval of 



b E (-V2a + 2, V2 a + 2) . 

At N > 2 a graphical determination of the domains D^ N > may be used. For 
illustration let us consider N = 3 and metric 

2 a 2 2«ia 

Q\'(a 1 ,a)= 2ga a + 1 ai 

with the «i— dependence of its three eigenvalues illustrated by Figure [1] at 
a = 1. 

The pattern of the graphical localization of the eigenvalues of our tridiag- 
onal metrics Q^\ai,a) remains qualitatively very similar in a broad range 
of parameters N, ol\ and a. In particular, we may be sure that the matrix 
®[ (ati,a) remains positively definite at all the sufficiently small nondiago- 
nalities, i.e., in a nonempty subdomain of D^ N ' where |«i| a. 

Several interesting as well as practically highly relevant questions arise 
when one tries to extend the graphical analysis to higher dimensions N. First 
of all, the growth of the necessary numerical precision makes the analysis a 
bit costly. Indeed, one must be careful with the numerical localization of 
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Figure 1: Three eigenvalues p = p(g) of metric Q\ (g, 1). 

the eigenvalues of the metric because even our fully explicit formula (j!7p 
leads to a perceivable numerical contrast between the maximal eigenvalues 
#o = 9\ — 1 and the unexpectedly quickly decreasing roots 9 7 ~ 0.0003968 
or 9 8 ~ 0.0000496 (etc) of the corresponding secular equation. 

Fortunately, the extremely elementary form of the matrix elements of our 
Gegenbauerian tridiagonal metrics (12 ip still supports the practical feasibility 
of the direct numerical localization of the boundaries of the related two- 
dimensional domains D[ N ^ of admissible parameters oc\ and a up to the fairly 
large dimensions. Moreover, there exists an encouraging numerical evidence 
that these boundaries dD^ stabilize and remain only very weakly dependent 
on the dimension at large N ^> 1. 

A persuasive sample of such an evidence is provided by Table 1 where 
we choose a = 1 and tabulated the values of G (our metric Q{ (g, 1) is 
positive definite for g G (—G,G)) together with auxiliary values G' and G" 
(our metric Q^\g, 1) has at most one or at most two negative eigenvalues 
in the larger intervals g G (— G', G') and g G (-G", G"), respectively). 

The Table strongly and very persuasively supports the N— independence 
of G ~ 0.776 in the limit N — > oo (i.e., the existence and stability of a non- 
empty domain d[ n ^ where the metric is positive). Indeed, the left column of 
the Table indicates that the N— th value of G only differs from its predecessor 
in the (N — 3)— rd decimal digit. 
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Table 1: The N— dependence of boundaries ±G of the domain D\ ' at a — 1. 



N 


boundary value G 


neighboring G' 


next G" 


1 


oo 






2 


1 






3 


0.8164965809 


— 


— 


4 


0.7835809235 


2.210430034 




5 


0.7772152453 


1.528761895 




6 


0.7761738933 


1.347821298 


3.702152325 


7 


0.7760367842 


1.284679682 


2.333798009 


8 


0.7760220038 


1.261982266 


1.922171587 


9 


0.7760206592 


1.254396565 


1.747726425 



6.3 Pentadiagonal metrics 

At k — 2 and variable N we may try to solve Eq. (|12p by pentadiagonal 
ansatz 
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(25) 



using the same recurrent method as above. The selection of 71 = a and 
the consistent specification of 72 = (1 + a)/(4 + 2a) initiate now the com- 
bined recurrences for two unknown sequences in f l25|) . These recurrences 
may be extracted, from linear algebraic Eq. ffl2|) . as a subset of all of its lin- 
early independent items. The result of their solution (which was, naturally, 
computer-assisted and rather lengthy) can be written down in closed form, 
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with j = 3, 4, . . . , N - 1 in 

?i = (o-lt\r ' T n = (1 + 2a) (2 + 2a) . . . (n + 2a) (26) 
(2j + 2a)r i _ 2 

and with To = 1 and j — 1, 2, . . . , N — 1 in 

c = 2[2a 3 + 3a 2 - (4j - 5)j a - (2j 2 - l)(j - 1)] 
J (2j + 2 + 2a)(2j-2 + 2a)r j _ 1 ' 1 j 

At the smallest subscripts j there occur incidental factorizations which sim- 
plify slightly the numerators, 

_ 2 (2a 3 + 3a 2 + a) _ (a + 1) (2 a + 1) 

1 ~ 2 a (4 + 2a) ~ 2 (a + 2) 

and 

2 (2a 3 + 3a 2 - 6 a - 7) 2 a 2 + a - 7 



(6 + 2a)(2 + 2a)(l + 2a) 2 (2 a + 1) (a + 3) ' 
The last missing element cu^ = u ( - N \a) in formula ( l25l) is exceptional. Due 

to its manifest truncation-dependence, its value must be computed, at each 

N > 3, by direct insertion in Eq. f lT2|) . At the first few dimensions this is an 

easy calculation which gives the (incidentally, negative though comparatively 

simple) series of formulae 

w (3) = - = 5 + 3q (28) 

2(1 + 2a)' " 2 (1 + 2a) (1 + a) (2 + a) "'" 1 ; 

Their extrapolation inspires the general ansatz 

( jv) = Kr+^jya) f . 

(2AA-4 + 2a)rjv_ 2 v ; 

and its subsequent confirmation giving 

u N = (2A-3)(A-2), v N = 3N-6. (30) 



This completes our closed-form construction of pentadiagonal solutions ([25 
of Dieudonne's Eq. (j!2p at any matrix dimension TV = 3, 4, 



7 Discussion 



In contrast to the recent theoretical experiments with discrete models pos- 
sessing point-like impurities [6] or boundary terms [13] , the interaction in our 
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one-parametric solvable toy model is a smooth function of position. This is 
an innovation which may be considered natural. In various limits and dy- 
namical regimes we may then specify energies E n and wave functions \ip n ) 
using the broad menu of formulae available for orthogonal polynomials in 
question. In our paper the eigenstates of H were selected, for the sake of 
definiteness, in the closed form of Gegenbauer polynomials. 

In the context of mathematics the main obstacle of calling the related 
solvable matrices H Hamiltonians appeared in their asymmetry (i.e., non- 
Hermiticity). This seemed to disqualify these matrices from playing the 
role of operators of observables. Fortunately, such a conclusion would be 
erroneous. The clarification of the paradox dates back to Scholtz et al [18] 
and Bender et al [32]. We just recalled and used their argumentation in a 
new concrete application. 

Our method of the reconstruction of the metric based on the use of 
discrete Hamiltonians and mediated by the computer-assisted solution of 
Dieudonne's Eq. f[T2l proved very efficient. It led to compact analytic for- 
mulae for a family of metrics. New discrete-lattice quantum model has been 
found as described by the pair of matrices (H, Q). The first component of this 
pair is the N— dimensional Gegenbauerian Hamiltonian H which has been 
chosen tridiagonal. The second component G of this pair is the reconstructed 
(and non-unique) metric. 

In a historical detour let us remind the readers that nuclear physicists 
opened this Pandora's box of Q 7^ / cca twenty years ago [IB] when consid- 
ering fermionic Hamiltonians f) (acting in complicated Fock's space of 
"physical" states \ip^)) as transformed into isospectral operators H (acting 
in another "friendly" space n {F) ). The net gain was that the bound-state 
energies became obtainable by the diagonalization of the simplified bosonic 
Hamiltonian H 7^ . The price to be paid was that the latter operator 
proved manifestly non-Hermitian in the usual, "friendly" Hilbert space 
with trivial G^ = /. 
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In other branches of physics the recipe has been revitalized in connection 
with the emergence of VT— symmetric quantum systems [55| I5H 155] . This 
opened new horizons in particle physics [5B] and in relativistic quantum field 
theory [57]. The key theoretical idea of the formalism (viz, the nontriviality 
of the product Q^Q :— O ^ I) remained the same but the philosophy has 
been changed. In place of starting from the knowledge of the physical, self- 
adjoint h = and from the subsequent clever choice of a simplifying map 
Q, the updated model-building strategy (cf. [28],|58]) takes a manifestly non- 
Hermitian "friendly" candidate for the Hamiltonian H = ^ (i?^) 

and tries to reconstruct the "physical" Hamiltonian f) = via Eq. ([TU]) . 

Our present proposal of a new solvable model was inspired by the main 
weakness of the latter scenario which lies in a huge uncertainty and ambi- 
guity of the assignment H — > f) marked, say, by an iV— component multi- 
index A attached to Q = f2(A)). This ambiguity was inessential during the 
nuclear-physics mappings f) — > H(X). In VT— symmetric context and in its 
pseudo-Hermitian generalizations [39J it is more serious. It implies the non- 
uniqueness of physics represented by the A— dependent operator f) ( A) . The 
same initial operator H admits many experimentally non-equivalent physi- 
cal interpretations. The variations of A generate non-equivalent self-adjoint 
Hamiltonians I) (A). This means that the same spectrum of energies may 
coexist with different observable characteristics (e.g., coordinates [5U| BP]). 

The suppression of the ambiguity of the multi-indexed mappings f2(A) 
and of Hamiltonians f)(A) may be performed, according to Scholtz et al. [18J, 
via an explicit specification of some other observables C, T> etc. They have to 
obey the same Dieudonne's conditions of cryptohermiticity. In practice, this 
goal may be achieved by requiring that one of the observables used for this 
purpose is a charge with involutivity property C 2 = I [2S]- In our present 
considerations we used another strategy proposed in Ref. [30] and based on 
the hypothesis of existence of a nontrivial, fundamental "smearing" length. 

We showed in J5U] that the smearing length does not vanish and does not 
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diverge in models where some of the the metrics possess the (2/c+l)— diagonal 
band-matrix form = 0^. The subscript k — 0, 1, . . . has been interpreted 
there as the measure of the size of the smearing. 

The simplest physical scenario of this form certainly emerges when one 
decides to use just the diagonal metrics ©o ^ I- In Ref. [30] as well as in 
our present concrete model this "no-smearing" option proved allowed. The 
related diagonal-matrix operator of the coordinate remained merely scaling- 
non-invariant. Our quantum Hamiltonians then became tractable as living 
on deformed but still local one-dimensional discrete N— site lattices. 

Once we turn attention to our present model and to its generic band- 
matrix metrics 0&(a) with 1 < k <C N, the picture is changed and the coor- 
dinates prove smeared [H]. This feature could make our elementary solvable 
model tractable, e.g., as a weakly and controllably non-local alternative to a 
deformed local k = lattice [501132]. 

On the experimental level one expects that such a weakly nonlocal sce- 
nario and its consequences (including, e.g., phase transitions) might find 
simulations in classical systems. A decisive theoretical as well as experimen- 
tal progress in this direction has already been reported in optics [I3JIII]. The 
practical implementation of the parallel experimenting in quantum world is 
hindered by several mutually interrelated obstacles. The most serious one 
may be identified with a certain conflict between the simplicity of the matrix 
H and the complicated guarantee of its Hermiticity via metric O. Our reso- 
lution of this conflict has been based on the simultaneous simplicity of both 
the operators H and O. 

The main theoretical profit provided by the fully non-numerical tractabil- 
ity of our model may be seen in its manifest compatibility with postulates 
of Quantum Mechanics in which one works, simultaneously [2J, with a triplet 
of Hilbert-space representations %( p > F ' 5 ) of the quantum system in question. 
The Hermiticity status of operators depends on the space but they only stay 
non-Hermitian in the "naive" and "false" space tv- F ^. Thus, in our model, 
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the knowledge of the friendly input matrix H ^ W is complemented by the 
equally friendly nature of the ad hoc metric 6 = 6^ 7^ / and, ipso facto, of 
the reconstructed standard Hilbert space %^ s \ 

In practical terms our Gegenbauerian example exhibits several specific 
friendly features. First of all, it is nontrivial that our metrics are banded. 
This property only followed from the explicit solution of the Dieudonne's 
equation. Secondly, the matrix elements of the pseudo-metrics (i.e., of the 
sparse-matrix components Vj of the metrics) emerged as elementary func- 
tions of the free real parameter a. Last but not least, the matrix elements of 
the diagonal, tridiagonal and pentadiagonal metrics exhibited even an almost 
complete independence of the truncation N. 

All of these features of our Gegenbauerian model reconfirm the feasibility 
of our original intention of finding a new model-building recipe. Certainly, 
this (and similar) solvable models would guarantee a viability of fitting many 
measured (and not just equidistant) iV— plets of levels ^ x P enmental ) ^he 
suitable N— plets ^ heorettcal ) Q f ^he well known zeros of an appropriate (i.e., 
in our exemplification, Gegenbauer) classical orthogonal polynomial. 
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Appendix A. Long-range metrics with k = N— 1 



Equations ([25]) + ([28]) with iV 



of the metric 0), with maximal A; 



3 offer the simplest nontrivial example 
N — 1 in which some of the matrix 
elements become truncation-dependent. We found that this form of man- 
ifest N— dependence characterizes all the Gegenbauer metrics with k > 2. 
In this sense, the diagonal and tridiagonal metrics appear exceptional. In 
principle, one could hope that a similar exceptionality could characterize the 
antidiagonal-like metrics which were found in some other models [30] and 
which could tentatively be characterized by the triangularity property 



for j < k . 



(31) 



The failure of these expectations can already be detected at the next dimen- 
sion N = 4 because the explicit violation of antidiagonality already charac- 
terizes the heptadiagonal pseudometric 



( 4 ), 












q 2 +2a+l 
a+3 







a 2 +2a+l 
a+3 



3a+5 



a 



3a+5 



(a+3) (2 a+1) 




(32) 



(a+3)(2a+l) 

We see that the loss of the up-down symmetry is transferred from the Hamil- 
tonian H to the metric. Thus, one can only expect that at a given N, the 
most elementary longest-range component Pjv^_\ (a) of the Gegenbauer met- 
rics will possess the following triangular equal-color form 



v { n n !m 



o 














pn 

pia 

P21 

.■ .■ 



: p u y y 

ft3 ft .■' 

Pl2 p 2 2 p 32 

Pll p 2 l p 3 i 



(33) 
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The quick growth of complexity of the, presumably, closed but much less 
compact formulae for the matrix elements in ([3"3l may be illustrated for in- 
termediate N = 8 for which the maximal-range 15-diagonal (pseudo)metric 
matrix may be constructed by solving Eq. (fl2|) via ansatz (1331 . In the nor- 



malization where 



= a our calculations yielded the elements 

17V J 



2a 2 + a (a + l)(a + 3) 2 a 3 + 9 a 2 + 10 a + 3 



Pn = a = — — — , pi2 



Pl3 



2a + 1 ' ^ a + 7 (2a + l)(a + 7) ' 

2a 4 + 17a 3 + 52a 2 + 67a + 30 _ (a + 1) (a + 3) (2 a + 5) (a + 2) 
(2a + l)(a + 7) (a + 6) " (2a + l)(a + 7) (a + 6) 

2 a 5 + 25 a 4 + 124 a 3 + 305 a 2 + 372 a + 180 



Pl4 



(2a + l)(a + 7) (a + 6) (a + 5) 
(2 a + 5) (a + 3) 2 (a + 2) 2 



(2 a + 1) (a + 7) (a + 6) (a + 5) ' 
plus perceivably less compact 

3 a 2 + 19 a + 26 (3 a + 13) (a + 2) 

P21 — — <J 77; , 1W — ; „n / — —7^: — _ <J 



(2 a + 1) (a + 7) (a + 6) (2 a + 1) (a + 7) (a + 6) 

3 a 5 + 52 a 4 + 342 a 3 + 1064 a 2 + 1551 a + 828 

P22 = — 5 ^ = 

(2 a + 1) (a + 6) (a + 5) (a + 1) (a + 7) 2 



(a + 4) (3 a 2 + 22 a + 23) (a + 3) 2 

-5 



\2 ' 



(2 a + 1) (a + 6) (a + 5) (a + 1) (a + 7) z 

6 a 7 + 143 a 6 + ... + 40218 a 2 + 37901 a + 14640 

P23 = —6 7. = 

(a + 5) (2 a + 3) (2 a + 1) (a + 1) (a + 7) 2 (a + 6) 2 

(2 a + 5) (a + 3) (3 a 5 + 55 a 4 + 380 a 3 + 1223 a 2 + 1811 a + 976) 
(a + 5) (2 a + 3) (2 a + 1) (a + 1) (a + 7) 2 (a + 6) 2 
(where the higher-degree polynomials in numerators have non-integer roots 

which are all real) and 

(a + 4) (3 a 5 + 25 a 4 - 78 a 3 - ... - 2025) 

P31 — — 2 3 ^ , 

(a + 5) (2 a + 3) (2 a + 1) (a + 1) (a + 7) 2 (a + 6) 2 

(where the fifth-degree polynomial in numerator has solely three real non- 
integer roots), 

[a + 3) (6 a 7 + 95 a 6 + 176 a 5 - 5106 a 4 - ... - 82280) 



P32 



2 (a + 5) (2 a + 3) (2 a + 1) (a + 2) (a + 1) (a + 6) 2 (a + 7)' 



2S 



(where the seventh-degree polynomial in numerator has solely five real non- 
integer roots) and, finally, 

6 a 8 + 59 a 7 - 621 a 6 - ... - 60120 

P41 — o 5" 

2 (2a + 5)(2a + 3)(2a + l)(a + 5)(a + 2)(a + l)(a + 6) 2 (a + 7) 3 

(with just four real and four complex roots of the eighth-degree polynomial 
in the numerator). Summarizing, these results demonstrate not only the 
efficiency of our computer-assisted algorithms but also, in parallel, the quick 
decrease of the practical appeal of working with more-than-pentadiagonal 
metrics Q^\a) with k ^> 2. 
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